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Simple wavelet sets in W 1 

Kathy D. Merrill 

Abstract. Wavelet sets that are finite unions of convex sets are constructed in R™, 
n > 2, for dilation by any expansive matrix that has a power equal to a scalar times 
the identity and also has all singular values greater than ^[n. In particular, we produce 
simple wavelet sets in any dimension for dilation by any real scalar greater than 1. 



< 

(~| 1. Introduction 

-i— > 

A wavelet set relative to dilation by an expansive (all eigenvalues greater than 1 

in absolute value) real n x n matrix A is a set W C lR n whose characteristic func- 

tion lw is the Fourier transform of an orthonormal wavelet. That is, if ip — lw, then 

[V>j,fc = Vl det A\Ji/j(A j ■ -k), J GZ.fce Z n | is an orthonormal basis for L 2 (R n ). This 

definition is equivalent to the requirement that the set W C W 1 tiles n— dimensional space 

(almost everywhere) both under translation by Z n and under dilation by the transpose 
CO 

^h A* } so that 

y lw(x + k) = 1 a.e. x G M n , and 

fcez n 

J^ l w (A* J x) = 1 a.e. x G M n . 

While wavelet set wavelets are not well-localized, and thus not directly useful for 
applications, they have proven to be an essential tool in developing wavelet theory. In 
particular, wavelet set examples established that not all wavelets have an associated MRA 
[S] , and that single wavelets exist for an arbitrary expansive matrix in any dimension [S] . 
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Smoothing and interpolation techniques have also used wavelet set wavelets to produce 
more well-localized examples. (See e.g. p3], [32], [7j , [6], p, [33], [3].) 

All of the early examples of wavelet sets for dilation by non-determinant 2 matrices 
in dimension greater than 1 were geometrically complicated, showing the fingerprints of 
the infinite iterated process used to construct them. (See e.g. Figure [l]( a)). Many early 
researchers, e.g. [6], |16| . conjectured that a wavelet set for dilation by 2 in dimension 
greater than 1 could not be written as a finite union of convex sets. In support of this 
conjecture, Benedetto and Sumetkijakan [6] showed that a wavelet set for dilation by 2 
in R n cannot be the union of n or fewer convex sets. However, in 2004, Gabardo and Yu 
|10] used self-affine tiles to produce a wavelet set for dilation by 2 in IR 2 that is a finite 
union of polygons (Figure [lib)). In 2008 [14J we used a technique based on generalized 
multiresolution analyses [2] to construct such wavelet sets for arbitrary real (d > 1) scalar 
dilations in M 2 . Figure [nc) shows one of the wavelet sets for dilation by 2 from |14j . 
Although they were developed independently, and using very different techniques, these 
two examples are remarkably similar. In fact, the wavelet sets in Figure [lib) andHTc) are 
equivalent in the sense that one can be transformed into the other under multiplication 
by a determinant 1 integer matrix. The similar shape of these two wavelet sets suggests 
the general n-dimensional result produced in this paper. 







^^^ (1.2, 



^^^^ . -1 , - . 4 
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Figure 1. Wavelet sets for dilation by 2 in M 2 

We call wavelet sets that are finite unions of convex sets simple wavelet sets. In 2012 
|15| . we expanded the results in |14j to produce simple wavelet sets for dilation by any 
2x2 matrix which has a positive integral power equal to a scalar times the identity, as 
long as its singular values are all greater than v2. In that paper, we also found examples 
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of expansive 2x2 matrices that cannot have simple wavelet sets. It is our conjecture that 
an expansive matrix whose determinant does not have absolute value equal to 2 can have 
a simple wavelet set if and only if it has a positive integer power equal to a scalar times 
the identity. 

In this paper, we generalize the 2-dimensional examples in [15] to n-dimensional space, 
n > 2. We do this using neither the generalized multi-resolution analysis techniques 
of |15j . nor the self-affine techniques of [10] . Rather, we use a remarkable result by 
Sherman Stein |17j on tiling IR n by notched cubes, together with the tiling conditions that 
are equivalent to the definition of a wavelet set. Section 2 skews and translates Stein's 
notched n-cubes to produce notched parallelotopes that are wavelet sets for dilation by 
negative scalars. Section 3 further modifies these notched parallelotopes by translating 
out a central parallelotope (as in Figure [l](b) and[I|c)), thus creating wavelet sets for 
positive scalar dilations, and more generally for matrices that have a positive integral 
power equal to a scalar, as long as their singular values are not too small. 

2. Tiling with notched parallelotopes 

We begin by establishing some notation. Write {ex, e^, • • • e n } for the standard basis 
of lR n , and C for the cyclic permutation matrix with columns (e2,e^, . . . ,e n ,ei). Let 1 
stand for the vector (1, 1, . . . , 1) e R n , and write r* for translation by t e M. n . 

Given a vector v in M. n that is not a multiple of 1, let 

V[v] = {x v + x x C(v) + ... x n - X C n - x v : < Xi < 1} 

be the parallelotope spanned by the vectors {v,C(v), . . . ,C n ~ 1 (v)}. Note that V[v] has 
two vertices on the line determined by 1: one at the origin, and the other determined by 
the sum of the coordinates of v , at (^ u»)l. Given an a, < a < 1, write J\f[v , a] for the 
notched parallelotope that results from deleting a subparallelotope scaled by a from the 
vertex (Y2 v i)l °f ^M- That is, let 

tf[v,a]=V[v]\ r (1 _ a)(5>)r ctf>M. 

We will need the following result about translation tilings by notched cubes due to 
Sherman Stein. 
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Lemma 2.1. Given a real number < a < 1, let L be the lattice spanned by the 
columns of I — aC, where C is the cyclic permutation matrix. Then the translates of the 
notched unit cube M[e\, a] by the vectors in L tile W 1 . 

Proof. See [T7] □ 

We use this result to produce a notched parallelotope that tiles W 1 under translation 
by the lattice Z n : 

Lemma 2.2. For a fixed real number a, < a < 1, letw(a) = x _ w (l, a, a 2 , . . . ,a n_1 ). 
Then the translates of J\f[w(a), a] by Z n tile M. n . 



Proof. By Lemma 2.1 we know that A/"[ei,a] tiles M. n under translation by L, the 
lattice spanned by the columns of I — aC. If we define A to be the linear transformation 
that maps L to Z n , we thus have that that A{J\T[e.\, a]) tiles M n under translation by Z n . 
Note that A = (I - aC)' 1 = J2Zo( aC Y = TZ^T^iotC) 1 . Thus, applying A to the 
notched unit cube Af[e\, a] will result in the notched parallelotope J\f[w(a), a}. D 

To be a wavelet set, a notched parallelotope would have to also tile under dilation; that 
is, its inverse dilate would need to fit perfectly into its notch. As this is clearly impossible 
with the origin as one of the extreme points of the parallelotope, we consider translates 
r t M[w(a),a\. For the dilate of T t M[w(a),a} by k, d > 1 to fit into its notch would 
require a = k and kt — t + 1, so that t = —-^1. Since this again forces the origin to be 
an extreme point of the parallelotope, tiling by dilation cannot work for positive scalar 
dilations. However, for negative scalar dilations, a wavelet set that is just a notched 
parallelotope is possible, as the following theorem shows. The examples produced by 



Theorem 2.3 generalize the wavelet sets for negative scalar dilations in IR 2 found in |10j 
and [T5] . 

Theorem 2.3. For d e R, d > I, let w{\) = ^i{d n , cT" 1 , ...,d), and t = —^. 
Then 



•ti 



dl'd 



is a wavelet set for dilation by —d in W 1 . 



PROOF. We know from Lemma 
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that N [w(i), 2] , and thus W; tiles M " under 



2.2 



translation by Z n . It remains to show that W tiles under dilation by — d. 

The parallelotope ^ (r t f P[w]) has its vertices on the line determined by 1 at ( j^) 1 
and at (-^zj) 1- Since the inside corner of the notch of r t j A/" [iu(i), |] is at (t + 1)1 = 
— :«it1> this shows that W U — ^W = r t j "P[iu] \ 4- (r t f -P[w]), which tiles under dilation 
by d 2 . Thus W tiles under dilation by —d. □ 



Figure^ shows the wavelet set produced by Theorem 2.3 for dilation by —2 in 



Figure 2. Wavelet set for dilation by -2 in 



3. Notched parallelotopes with satellites 

Even though we have seen that a translate of a notched parallelotope r< J\f[w(a),a] 
cannot itself tile under dilation by a positive scalar, we will make a simple alteration to 
such a set that retains the property of tiling under translation, and makes the set tile 
under dilation as well. We use the idea behind the wavelet set construction technique in 
j3]. That is, we eliminate overlap between our proposed wavelet set and its dilate, by 
translating the dilate out by an integer vector, creating a satellite. In order to avoid the 
iterated process required in [4], we translate out a little bigger piece than the dilate of 
Tt Af[w(a), a}; that is, we translate out the dilate of the whole parallelotope r t V[w]. We 
choose the translation amounts such that a dilate of the satellite exactly fills the notch. 
The details of this construction are carried out in the following theorem. 
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Theorem 3.1. Forde R, d> 2, letw{^) = ^^{d 2 ™ , d 2n ~ 2 , . . . ,d 2 ). Suppose k e Z 
satisfies 1 < k < d, and let t = jp", • Then 



W = I T tt N 



ir[ dV^ 



\ (^fW) \Jr kf (^r tT V[w} 



is a wavelet set for dilation by d in 



Proof. We claim that \r t i V[w] C r t y A/"[w(^),^]. To see this, first note that 
t < and — £ < -^— j- = X^=i w i> so that e ^ r ti ^M C r t j P[w]. Thus, to establish 
the claim, it will suffice to show that the vertex of the notch that is closest to the origin, 
namely (t + (1 — jp)Yli w i)^-i ^ es outside of 4r t f V[w}. That is, we must show that 
t + 1 > i It + ^— j- j. Substituting t = ^ 2 ~ , we see that this is equivalent to fc > ^-j-, 
which follows from the given conditions k > 1 and d > 2. 

impies that A/" [w(^), ^] tiles under translation by the integer lattice, and 



Lemma 



2.2 



thus that r t f M [w(^), jz] does as well. By the claim, we have that W is formed from 
T t fN ju>(Jj), Jjl by translating the subset hr t f ^M by an element of the integer lattice. 
Thus W tiles lR n under translation by Z n . 

To establish tiling under dilation by d, we first show that T t fN [^(Ja), ji] is disjoint 
from r k i (h^ti ^M)- That is, we must show that 4 + k > t + 1, which follows from the 
definition of t together with the condition k > 1. Now, note that r t j T^to] \ -jT t f V[w] tiles 
1R™ under dilation by d. The definition of t together with the disjointness of the pieces 
of W shows that W is formed from this set by dilating the notch T, t+1 ^ (j-z'Plw]) by d. 
Thus, W also tiles lR n under dilation by d. D 



Remark 3.2. The wavelet sets produced by Theorem 3.1 are a natural generalization 
of the 2 dimensional simple wavelet sets produced for scalar dilations in |14] . We can also 
alter these examples to produce natural generalizations of the 2-dimensional example for 
dilation by 2 that appears in [10J. Note that if W is a wavelet set for dilation by the 
scalar d, then so is S(W) for any integer matrix S of determinant ±1. If we take S to be 
the n x n matrix that has l's on the diagonal, -l's on the subdiagonal and O's elsewhere, 
then S(W) is a simple wavelet set for dilation by d in dimension n that is centered on the 
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X\ axis rather than the line x\ — x-i — . . . x n . Other variations are easily produced using 
other choices for the matrix S. 



Figure [3] below shows one of the wavelet sets produced by Theorem 3A for dilation by 



2 in 



as well as a variation described in Remark 3.2 





Figure 3. 3-dimensional wavelet sets for dilation by 2. 



The final theorem uses the same technique as Theorem 3.1 to produce simple wavelet 



sets for dilation by expansive matrices that have a power equal to a scalar multiple of the 
identity, as long as their singular values are all larger than \/n. 

Theorem 3.3. Let A be an n x n integer matrix such that A p = d(I[n}), where d > 1, 
p G 7L, p > 1, and I[n] is the n x n identity matrix. Suppose further that all the singular 
values of A are greater than y/n. Let w(-h) = dq ^_ 1 (d qn , d qn ~ q , . . . , d q ). Let k be the closest 
vector in Z n to A*~ l ( yl ) , and let t = -^—^ ( A*k — d q l ) . Then for q sufficiently large, 

-) -' 

d« r d q 



w 



n 



M 



w 



\ A*- l r t V 



XV 



s)])u(^(^"Ks 



is a wavelet set for dilation by A in 



Proof. By Lemma 



2.2 



M[ 



w{ 



l_ 

dt)i di 



4-1 , tiles under translation by Z™, and thus its 
translate by t does as well. Thus, W will also tile by translation as long as the set 
A*~ l r t V[w(-^)], which is translated out by k, is a subset of T t M [w(-^), j^]. As q — > oo, 
the set t{P [w(^-)1 approaches the unit n-cube centered at the origin. Thus, by taking q 
sufficiently large, the longest vector in r t V [w(^)] will be arbitrarily close to y/n times 
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as long as the shortest vector. Then, since the singular values of A are greater than s/ri, 
we will have A*^ 1 [r{P [w(J?)]) C r{P [w(-4-)l . The size of the notch also shrinks to as 
q — > oo, so that by taking q larger if necessary, we will also have A*' 1 (r t V [w(^)]) C 

To establish tiling under dilation, note that because of the containment established 
above, r t V[w{±)\ \ (A*- l r t V[w(±)\) tiles W l under dilation by A*. The definition of 
t assures that the outlier piece of W exactly fits into the notch of r t H [w(-h), 4-1 under 
dilation by A*^ pq+1 \ Thus we have that W also tiles M. n under dilation by A*. □ 



Remark 3.4. Theorem 3.3 also produces simple wavelet set for scalar dilations 1 < 



d < 2, which were not covered by Theorem 3.1 For scalar dilations by d > 2, Theorem 



3. 3| produces a series of alternative wavelet sets to those of Theorem 3.1 As q increases 
in this series, the parallelotope becomes closer to cubic, the notch becomes smaller, and 
the satellite becomes farther removed. 



/ 



Example 3.5. Let A* 



3 

3 

1 -3 



\ 



/ 



Then A 2 = 91 and SingularValues^) 



{3.54, 3, 2.54}, so Theorem 3.3 applies. With w(|) 
and t = (-|,-f,-§), we have A^nV [w(|)] C r t V [ 
wavelet set 



9- 



!') 



W 



9 3 -l' 9 3 -l' 9 3 -l/' 

1 



k= (1,1,-1) 



. Thus, we have a simple 



W 



T t 



M 



w 



-) - 

9 '9 



\ A*-\ V 



IV 



G)])u(^w,(^^[-g) 



which is pictured in Figure 4. 



The hypothesis in Theorem 3.3 that the singular values of A be greater than y/n is 



sufficient but not necessary, as the next example shows. 



Example 3.6. Let B* 



( 



2 





1 





-2 











-2 



{2.56, 2, 1.56}, so Theorem 3.3 does not apply. With w{\ 

i 



i\ 



/ 



. Then B 2 = 41 and Singular Values (B) 



and t = (— 1, — |, — |), we do not have B* 1 r t V [ 



w( 



i 

C T, 



'64 16 M L 
-63' 63' 63 J' 



v[ 



w{ 



;i, -l, -i) 

(See Figure 5(a).) 
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Figure 4. Simple wavelet set for the matrix A of Example 3.5 



However, using q = 2, with tu(i) = (^ri, ^nnri' ^=l)> * = ~ll^ and fc = ( 6 ' ~ 4 ' ~ 4 )' 
the required containment does hold, yielding the wavelet set 



W 



n 



N 



w 



1 

16 i ' 16 



\ B-V t P 



■u> 



16 



|J fr k (b*-\ V 



w 



16 



Figure 5(b) shows the central part of the wavelet set. (The complete wavelet set also 
includes a translation of the missing inner parallelotope by (6, —4, —4).) 





(a) B- l V [w (I)] t V [w (J)] 



(b) ^> (£),&] \ B- 1 ? K£)] 



Figure 5. Building a simple wavelet set for the matrix B of Example 3.6 
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